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Many problems can be solved by preparing a specific eigenstate of some Hamiltonian H. The 
generic cost of quantum algorithms for these problems is determined by the inverse spectral gap of 
H for that eigenstate and the cost of evolving with H for some fixed time. The goal of spectral gap 
amplification is to construct a Hamiltonian H' with the same eigenstate as H but a bigger spectral 
gap, requiring that constant-time evolutions with H' and H are implemented with nearly the same 
cost. We show that a quadratic spectral gap amplification is possible when H satisfies a frustration- 
's^ free property and give H' for these cases. This results in quantum speedups for optimization 
problems. It also yields improved constructions for adiabatic simulations of quantum circuits and 
for the preparation of projected entangled pair states (PEPS), which play an important role in 
quantum many-body physics. Defining a suitable black-box model, we establish that the quadratic 
, ^ amplification is optimal for frustration-free Hamiltonians and that no spectral gap amplification is 
possible, in general, if the frustration- free property is removed. A corollary is that finding a similarity 
transformation between a stoquastic Hamiltonian and the corresponding stochastic matrix is hard 
in the black-box model, setting limits to the power of some classical methods that simulate quantum 
adiabatic evolutions. 
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^ I. INTRODUCTION AND SUMMARY OF RESULTS 

^ Many problems in physics and optimization, such as describing quantum phases of matter or solving satisfiability, 

qh can be reduced to the computation of low-energy states of Hamiltonians (see, e.g. [121 [31]). Methods to compute such 
I— I states exist, with adiabatic state transformation (AST) being one of the most acknowledged and powerful heuristics 
for that goal [151 - 501 . The problem in AST involves transforming a pure quantum state |'(/'o) into these 

being the endpoints of a state path \^ps)^ < s < 1, where each \^s) is an eigenstate of a Hamiltonian H{s). Quantum 
strategies to solve the AST problem typically require evolving with the Hamiltonians for time T > 0, which determines 
the cost of the method. A well-known example is given by the quantum adiabatic theorem and requires changing the 
interaction parameter s sufficiently slowly in time [531 [2Z] j but more efficient strategies exist [HI [TOl [13] ■ When direct 
access to Hamiltonian evolutions is not possible, such evolutions can be simulated in the standard quantum circuit 
model by means of Trottcr-Suzuki-like formulas [HlUHlSll- For these cases, the cost of the method does not only 
y—^ depend on T, but also depends on the additional cost of approximating fixed-time Hamiltonian evolutions by quantum 
T-H circuits. For this reason, it is standard to restrict to those H{s) whose fixed-time evolutions can be approximated by 

quantum circuits of small size (e.g., size bounded by constant-degree polynomials). 
^ Let L be the angular length determined by {|V's)}o<s<i and A^. > A the spectral gap to the other nearest eigenvalue 

of H{s). Recently, we provided fast quantum methods for eigenpath traversal that do not exploit the structure of the 
Hamiltonians and solve the AST problem in time T G r2(L/A) [TD]. Such methods were shown to be optimal in a 
black-box model [TT. However, the results in Ref. [TT] do not forbid the existence of faster methods for solving the 
AST problem if additional knowledge about the structure of H{s) is available. Our paper focuses on this observation 
and introduces the spectral gap amplification problem or GAP. Basically, we intend to answer: Can we construct a 
Hamiltonian H'{s) that also has \ips) as eigenstate, but spectral gap A'^, ^ A^? 

When spectral gap amplification is possible, it can be used as a technique to find new quantum speedups in both, 
the Hamiltonian-based model and the quantum circuit model. In Refs. |341 136j. for example, we constructed the so- 
called quantum simulated annealing algorithm that provided speedups of the well-known simulated annealing method 
implemented using Monte-Carlo. The reason for such speedups is a gap amplification step: if A is the spectral gap of 
the stochastic matrix S used in Monte-Carlo, there is a Hamiltonian H' with spectral gap A' > \/A and eigenstate that 
allows to sample from the fixed point of S. To build H' we used previously-developed tools for quantum walks fS] I38j . 
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Motivated by the results in Refs. [Ml [33] we study the GAP in different scenarios. We first show that, for Hamil- 
tonians that satisfy a frustration-free property [131 Ull IMl HI] , a quadratic spectral gap amplification is possible and 
optimal in some suitable black-box model. Quadratic spectral gap amplification was previously known only for Hamil- 
tonians that result from a similarity transformation of stochastic matrices [3U[3ni[3H]- These Hamiltonians, which 
are the so-called discriminants, are also stoquastic, i.e. the off-diagonal entries are non-positive. A direct consequence 
of our new construction is that spectral gap amplification can now be achieved for any frustration-free Hamiltonian 
as those used in quantum adiabatic simulations of quantum circuits, improving upon the results in Refs. [Tl [?1 [2H]. 
The eigenstate in our construction is not the lowest-energy state of the Hamiltonian induced by the quantum circuit 
(which was the case in Refs. [T1[31I2H])- Nevertheless, techniques for AST work for any eigenstate and still apply to 
our case, giving more efficient quantum adiabatic simulations to prepare the circuit's output state. In addition, some 
low-energy eigenstates of general frustration-free Hamiltonians (termed PEPS) were shown to play an important role 
in physics, renormalization, and optimization [501 1301 HI] : PEPS arc useful variational states reproducing the local 
physics with high accuracy [33 HO] • Our method for spectral gap amplification can be used to speedup the preparation 
of PEPS on a quantum computer, with an expected improvement of other methods for this goal (e.g., [32]). 

Interestingly, it is hard to achieve spectral gap amplification for general Hamiltonians that do not satisfy the 
frustration-free property. It would otherwise imply fast quantum methods for some problems, contradicting known 
complexity bounds in the oracle setting |7 . In addition, a natural generalization of our technique for spectral gap 
amplification would output a Hermitian operator only if the input Hamiltonian H is semidefinite positive (as is the 
case for frustration-free Hamiltonians). Such a generalization could then be used to solve the decision problem of 
whether a given H is semidefinite positive or not (with some small error bound), and provide an estimate to the 
lowest eigenvalue of H. This is a complete problem for the class QMA [26], and finding its solution may require a 
quantum circuit of undesirably large size. 

The result on the impossibility of general spectral gap amplification has additional implications on the power 
of classical methods that solve the AST problem. For example, some stoquastic Hamiltonians can be mapped to 
stochastic matrices through a similarity transformation [12l [131 El] . The fixed point of the stochastic matrix coincides 
with the distribution induced by the ground state of the Hamiltonian. The stochastic matrix would allow us to build 
a Monte-Carlo method to sample from the fixed point and classically solve the AST problem in this case. The same 
similarity transformation could then be used to give a frustration- free representation of the Hamiltonian [331 140| . 
allowing us to use the results for frustration-free Hamiltonians to amplify the gap. This contradicts the impossibility 
result, implying that finding the transformation is hard (in the black-box setting). 

We organize the manuscript as follows. First, in Sec. [H] we define the GAP in more detail and present some methods 
for the proofs of the quadratic amplification in the frustration-free case, and for the impossibility of amplification in 



general. In Sec. Ill we show the construction of H' for frustration-free Hamiltonians. In Sec. HI A we prove that 



the cost of implementing evolutions under H' is nearly the same as the cost of implementing evolutions under H 
in the black-box model, and present a simulation method. In Sec. |IIIB| we prove that the quadratic amplification 



is optimal for frustration-free Hamiltonians in this black-box model. In Sec. IV we comment on the advantages of 
our gap amplification construction for the adiabatic simulation of quantum circuits and give a local Hamiltonian. In 
Sec. |V]we prove that no gap amplification is possible for general Hamiltonians in the black-box model. This results 



in limits on the power of classical simulations for quantum systems, that are discussed in Sec. VI We conclude in 
Sec.lvni 



II. DEFINITIONS AND METHODS 



We first comment on the implementation cost of Hamiltonian evolutions. Ignoring precision, we assume Hamil- 
tonians H given as black boxes that implement the evolution exp{—iHt} on ancilliary input \t), satisfying \t\ < tt. 
When we write H = Afc for the Hamiltonian, we assume a modified black-box that implements exp{— iAfci} on 
ancilliary input Such a black box is typically an efficient quantum circuit in applications. We write Oh for the 

black-box. The implementation cost of exp{— i_ff's} in the black-box model, where H' is the Hamiltonian resulting 
from the gap amplification construction applied to H, will be the number of times Oh is called to implement the 
operation in either case. 

Definition 1 Let H = X^fe^fe ^ ^ finite- dimensional Hamiltonian, \il>) a (unique) eigenstate of H with 

eigenvalue \, and A the gap to the other nearest eigenvalue of X; i.e. A is the spectral gap. The goal of the GAP is 
to find a new Hamiltonian H' that has \ip) €5 |o) as (unique) eigenstate and spectral gap A' G il{A^~^), for e > 0. |o) 
is any subsystem's state for which an efficient quantum circuit is known. The implementation cost o/exp{— in 
the black-box model must be of the same order as the implementation cost o/ exp{— i_fft}. 
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The last requirement forbids naive constructions such as H' = cH for c 1 and is needed to carry quantum speedups 
from the Hamiltonian-based model, in which the cost is the evolution time, to the black-box model, in which the 
cost is the number of black boxes. If both H and H' are sparse with a bounded number of (efficiently computable) 
non-zero entries per row and matrix oracles for the are provided, the last requirement can be satisfied pi 151 1151135) . 

Definition 2 A Hamiltonian H G x is frustration free if it can he written as H — ^fellfe , with (known) 

< flfe < 1, (Ilfe)^ = life projectors, and L G 0[polylog{N)]. Further, if [ip) is a lowest-eigenvalue eigenstate (ground 
state ) ofH, then Ilk\tp) = V fc. 

lip) is then a ground state of every term in the decomposition of H. We will assume that Ofe = 1 Vfc when there is 
no loss of generality. Then H = life is still frustration free, it has |?/') as ground state, it has at least the spectral 
gap of flfellfe, and ||ff || < L. While the life's may be local operators, we do not make that assumption here. It will 
suffice to have access to the black box Oh- 

We consider the unstructured search problem or SEARCH to prove the optimality of quadratic gap amplification for 
frustration-free Hamiltonians and also the impossibility of gap amplification for general Hamiltonians. We introduce 
SEARCH and describe a simple variant of Grover's algorithm to solve it. 

Definition 3 Consider a family of Boolean functions or oracles {fx}, with domain [0, ..,M — 1], that satisfy fx{y) — 
5xy The goal of SEARCH is to find the unknown word x, with high probability, by querying fx the least possible 
number of times. 

For quantum algorithms, we use a reversible version of fx defined as 

QAy)\b) = \v)\h®fx{y)), (i) 

with h G {0, 1}. We refer to Qx as the search oracle. Quantum algorithms that solve SEARCH with Q{^/M) uses of 
Qx, and other x-independent operations, exist [3 HI]- Classical algorithms require 0(M) queries to fx- 
From Ref. |16j we can give a quantum algorithm to solve SEARCH using projective measurements. 

Lemma 1 Let be a family of Hamiltonians with unique eigenstates \ipx) satisfying px = \ {x\ipx)\'^ G ^(1) ^"^"^ 

Ps — \{s\'tj)x)\'^ G ^(1); where \s) — -^jj'^y^ti \y) equal superposition state- A quantum algorithm can find x 

with cost 0{1/A), where A is a bound on the eigenvalue gap of Hx- 

Proof 1 The probabilities Px and Ps are bounded independently of M. If we initialize a quantum computer in \s) and 
next we perform two projective measurements of \il^x) o.nd \x) , respectively, we can find x with probability no smaller 
than Px-Ps- The cost of this method is measured by T, the time of evolving with Hx- This is needed to simulate 
the measurement of l^l^x), with high accuracy, using the well-known phase estimation algorithm \16l \25^ or by phase 
randomization J^. The cost T is dominated by the inverse gap of Hx- In optimal constructions, A G ©(l/V^M) so 
that T G 0{VM). 

III. GAP: FRUSTRATION-FREE HAMILTONIANS 

When a Hamiltonian H = J^k is frustration free, local, and stoquastic [T^, we could use the results in Ref. |T3] 
and Refs. [34l[36] to build a Hamiltonian H' with a quadratically bigger spectral gap. In that case, H is the so-called 
discriminant of a Markov process or stochastic matrix. Such a construction uses tools from quantum walks |38j to 
speed up Markov chain based algorithms and requires full knowledge of each Hfe. Here, we present a novel and 
more efficient technique that allows for quadratic gap amplification of any frustration-free Hamiltonian, without the 
requirement of H being local and stoquastic. Later we show that this technique provides the optimal amplification 
and present a black-box simulation method for the evolution with H' that does not require full knowledge of each Hfe . 

Tiieorem 1 Let H — X]fc=i ^fe G x be a frustration-free Hamiltonian, a ground state, and A the spectral 
gap or smallest nonzero eigenvalue. Then, there exists H' satisfying H'\ip)\o) = and eigenvalue gap A' G AL). 
The dimension of the null space of H' is that of the null space of H . 

Proof 2 The proof is constructive. We let 



L 

X = ^nfe®Tfe 
fe=i 



(2) 
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be a Hamiltonian, and = a rank-one projector acting on C^. It is easy to show that X" = X for n > 1. 

We define the unitary operator 

L 

(7 = e-'"^ = ]l-2^nfc0Tfc , (3) 

k=l 

with 1 the identity matrix whose dimension will be clear from context (here it is NL x NL). The implementation ofU 
requires a single black box Oh and = l. In addition, |o) = Sfe=i 1^) '■^ '^'^ ancilliary state and P = l (g) |o)(o| 
its projector. Because L € 0[polylog{N)], quantum circuits that prepare |o) efficiently exist (e.g. \35f). 
We define the Hamiltonian 

G = UPU -P . (4) 

If {\'^j)}i<j<N the eigenstates of H with eigenvalues Ai = < A2 < . . . < Ajv l£ L, G leaves invariant the at most 
two-dimensional subspace Vj spanned by {\ipj)\o) ,U\ipj)\o)}, for all j . To show this, we need an important property 
that is a consequence of our construction: 

PUP=(^l-^H^<S)\o){o\=A(»\o){o\. (5) 

The eigenvalues of A are jj — 1 ~ 2Xj/L and the eigenstates are also the To avoid technical issues we assume 

that "fN > 0. // this were not the case, the projector P could be easily modified so that 7j = 1 — Aj/L instead, and the 
assumption is fulfilled. We can think of as being cosaj = {ipj\{o\PU P\il>j)\o) , so that aj is the angular distance 
between the states \tpj)\o) and U\ilij)\o) (Fig. [7p. When j is such that Xj — 0, we have U\iljj)\o) = | "0^)10) and Vj 
is one dimensional. This implies that \^pj)\o) is an eigenstate of G of eigenvalue 0. When Xj ^ 0, Eq. ^ gives 
G'|'0j)|o) = jjU\'iJjj)\o) — \ipj)\o), and G is invariant in Vj. 



sinoj costtj -Lj) 










{UPU-P)\ib,)\o) 
















■* ■ 


/ ^2 





l^j>|o) 



FIG. 1: Geometric representation of the action of G on |i/;j)|o). 

When Xj 7^ (or OLj ^ we define |_Lj) as the unit state orthogonal to IV'j)!") in the subspace Vj. We can use 
simple geometry arguments (Fig.^ to show that the two-dimensional block of G in the basis {\^j)\o), |-Lj)} is 

„ f — sin^a, sin a," cos a 
Gj ^ [ . .■'2 

ysm aj cos aj sm Uj 

The eigenvalues of G in this subspace are ± sin aj . Because sin aj > y^2Xj/L, a quadratic amplification in the spectral 
gap follows. 

We need to account for a small technical issue to obtain H' . G acts trivially on any other eigenstate orthogonal 
to (BjVj. The dimension of the null space of G is much larger than the dimension of the null space of H . Usually, 
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this poses no difficulty in the AST problem because G does not generate transitions between the states \'il'j)\o) and 
any other state in the null space of G. In those cases we can assume H' = LG so that \\H'\\ < L. Otherwise, since 
P|(/i) — \''Pj){''Pj\ I'') ("I )!'/') = 0, we can easily remove the unwanted degeneracy of G by adding a "penalty term" 
to those states whose ancilliary register is not in |o). That is, 

H' = LG + 5{1 - P) (7) 

for some (5 > determined below. Note that, if Xj — 0, |i/'_,)|o) is still in the null space of H' . For Xj ^ 0, the matrix 
of H' in Vj is 

/ • 2 

^ ysin aj cos aj S/L + sin^ a. 

It suffices to choose S = \/ AL to guarantee that , the eigenvalues of H' in this subspace, satisfy |^*| > (1/6) -^/XJX. 

The gap of H' is then A' G 0(V AL) and the dimension of the null space of H' coincides with that of the null space 
ofH. 

Our choice of S uses A, which may be unknown. We can then safely use a lower bound on the spectral gap of 
H , A < A, and build H' using A instead. The loss is that the resulting gap of H' may not be the largest possible. 
However, this is not a drawback if we use H' (or LG) for the AST problem. As mentioned, the relevant spectral gap 
is of order V AL even if A is unknown, and the role of A is unimportant in this case. 

A. Black-box simulation 

In a Hamiltonian-based model of quantum computation, a gain on the spectral gap readily provides a quantum 
speedup for the AST problem. This model assumes the ability of evolving under H' , i.e. applying V{t) = exp{—ill't} 
directly. To carry those speedups to the quantum circuit model, we need a way of approximating a fixed-time evolution 
by a short sequence of gates. Usually, such a sequence is obtained by using Trotter-Suzuki- like approximations [511^. 
In those approximations, the evolution of a Hamiltonian H = A^ is built upon short-time evolutions under each 
Afc, and each such evolution is written as a constant-size sequence of gates |35) . In other words, what determines the 
size of the circuit that approximates V{t) is the implementation cost as given by the number of uses of Oh, the black 
box described in Sec. |llj Here we describe a simulation method for V{t) whose overall cost is almost linear in L\t\. 
In most cases this result implies a quantum speed-up, in the circuit model, for those methods that solve the AST 
problem using paths of frustration-free Hamiltonians. 

Lemma 2 For any k and e, there is a black-box simulation W(t) of V{t) that requires using the black box 
O times, k e N*, and satisfies \\W{t) - V{t)\\ < e. 

Proof 3 Up to a trivial offset, H' = Ai - {L + 6)P, with Ai = LUPU, \\Ai\\ = L, and \\P\\ = 1. We can then use 
known results on Hamiltonian simulations 1^5). \29^ to approximate exp{— iiJ't} at precision e by a sequence of 
c(e, K)|it|^^^/'^ concatenated evolutions under Ai and P. For fixed e and k, c(e, k) is constant. The proof follows by 
noticing that 

i.e. it requires two uses of U , and each U can be implemented using Oh once on input |7r). Evolutions under P for 
any time do not require the black box and can be implemented with a quantum circuit of size L or smaller. 

B. Optimal amplification 

We are interested in finding the biggest spectral gap amplification possible for frustration-free H. To this end, we 
consider the set T of all those H' that have I^A) |o) as eigenstate and such that evolving with H' for time t can be done 
using 0(1^1) black boxes Oh', i-e. L is constant. Under these assumptions, the quadratic spectral gap amplification is 
optimal in the black-box model: 




Theorem 2 A' G Q{VA). 
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Proof 4 We sketch the proof here and leave the full version in Appendix \a\ We consider instances {Hx}, x € 
{0, 1, . . . Af — 1} unknown, that can be used to solve SEARCH as described in Lemma^ Each "is frustration free 
and the spectral gap is of order 1/M . The corresponding black box Ofj can be implemented with a single call to the 

search oracle of Eq. If H'^ is the Hamiltonian with the amplified gap A', SEARCH can be solved by evolving with 
H'^ for time Q{l/A'). By assumption, this would require using fl(l/A') black boxes Oh^ and thus fl{l/A') search 
oracles. The lower bound on SEARCH implies A' e QiVA), with y/A e 0{1/\/M). 



IV. GAP AMPLIFICATION AND ADIABATIC QUANTUM COMPUTATION 



A property of our spectral gap amplification method is that it yields quantum speed-ups for the AST problem. It 
is also important to describe additional advantages of the method, including the simulation of quantum circuits. It 
is known that any quantum circuit specified by quantum gates Ui, . . . , Uq can be simulated, in a Hamiltonian-based 

model, by an adiabatic quantum evolution involving frustration-free Hamiltonians H{s) = X]fc=i '^fe('S)nfc(s). For 
s = 1, the ground state of H[l) has large overlap with the state output by the quantum circuit. In some constructions, 
L g C'[poly((5)] and Ilfc denotes nearest-neighbor (two-body) interactions between spins of corresponding many-body 
systems in one or two-dimensional lattices [U [3] . 

The cost of the adiabatic simulation depends on the inverse minimum gap of H{s), A. We can therefore use our 
gap amplification construction for frustration-free Hamiltonians to speed-up the adiabatic simulation. However, in 
doing so, we encounter difficulties regarding the locality of H' . For example, the projector P can be interpreted as 
an interaction term involving log2(i) qubits and this is undesired if we want to find a physical H' whose interactions 
address a limited number of subsystems. A simple procedure to avoid many-body terms is to consider each state \k) 
in |o) as a state of L qubits of the form 



\k)^\0...0 ^ 0...0) 

kth position 



(10) 



The subspace defined by Eq. (10) is typically termed the single-particle subspace. The terms appearing in P 

map to (Tj!"(Tj7,, with ct* = ((T^±icr^)/2, and a" the corresponding Pauli operator. The terms |fc)(fc| map to {l — a^)/2. 
In this subspace. 



G 



1 



fc(l-'^D/2^+^-p-"rnj,'(l-<')/2 



k,k' 



(11) 



If the Hk are two-subsystem interaction terms, H' — LG + (5(11 — P) will map to a combination of four-subsystem 
terms. This makes H' local but not spatially local. We also note that, while the overall spectrum of H' may have 
changed because of the mapping, the spectrum of H' in the single-particle subspace corresponds to the one analyzed 
in Sec. HIl 

We note that our construction to amplify the gap is not unique. Indeed, as seen in Appendix [Bl the Hamiltonian 



G = ^Hfe®(|fc)(0|-|-|0)(fc|) 



k = l 



(12) 



has a similar property: its eigenvalues are iy'Xj. Note that if a/c 7^ 1 we can modify it as G = Hfe 
(|fc)(0| + |0)(fc|). The action of G in the single-particle subspace is equivalent to that of j3S] 



fe=i 



(13) 



To eliminate the degeneracy of the null eigenvalue in this subspace, we can proceed as in Sec. Ill and add a term 
proportional to (1 - |0)(0|) ^ (1 + cr^)/2, so that 



G ^ G = ^ Hfe ® (a+ao + a,. C7+) 



+ , , %/A(l + ag) 



(14) 



k=l 



The relevant spectral gap of G in the single-particle subspace is still of order V A. 
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Nevertheless, an unwanted degeneracy of the null eigenvalue from other (many-particle) states is possible. (Note 
that any eigenstate in the 0-particle subspace has eigenvalue exactly •\/A/2). To remove such a degeneracy we set 
"penalties" to states that belong to different particle subspaces. To this end, we note that 



< G2 = ^ nkllk-{a+a^,a,a+ + H.c.) 
H.c. , 



fc.fe' 



< S S ctq (Tq 



(15) 



where — X]fe=i "^sS^ proportional to the so-called SU{2) spin raising operator. Let \S,m) be the eigenstates 
of = X]fc=i^fci with eigenvalue m, and total SU{2) spin proportional to S. Then, m e {— L,— L + 2,...,L}, 
< S < L, and 



S+\S,m-2) =N{S,m)\S,m) , 

NiS,m) = ^^{S + m)iS-m + 2) 



It implies 



{S+S-aoa+ + H.c.) |6)o ® \S, m) 



(16) 



(17) 



= {Sb^NiS, m + 2)2 + 6b^iN{S, mf) \b)a ®\S,m) , 



where b G {0, 1} denotes the state of the qubit in the 0-th position. The operator leaves the a-particle subspace 
invariant, and we write a = (L — m)/2. From Eqs. ( 16 1 and ( 17 1, the eigenvalues of G^ in that subspace are bounded 
from above by {L + m + 2){L — m-\-2) / A = (L -I- 1 — ajjl -I- a) . The upper bound to the absolute value of the eigenvalues 
of G, in that subspace, is ^J{L + 1 

WedefineZ = L-2-X;^ 
We also define 



a)(l + a). 

k=o '^k ^^'^ note that Z acts trivially in the single-particle space and nontrivially elsewhere. 



H' ^-G 
L 



2Z 



d > 0. H' acts invariantly in the a-particle subspaces. The eigenvalues of H' are then bounded from below by 



1 



^J{L + l-a){l + a) + 4(a - 1) 



(18) 



(19) 



For the specific case of d = 2, a = 1, we already know that the eigenvalues of H' are in the range [—1,1]; this 
case corresponds to the single-particle one whose spectrum was analyzed previously. Also, for d < 2, a > 2, all the 
eigenvalues coming from the a-particle subspaces do not mix with the eigenvalues of the single-particle subspace, and 
the potentially high degeneracy of the null space is avoided. For a = the exact eigenvalues of H' are \/A/ (2L^/'')— 4 < 
-7/2 if d < 2. 

Summarizing, the Hamiltonian 



H' = 



1 



^nfc(o-+cro 



A(l + a°) 



.k=l 



2Z 



(20) 



has |V'o)|10 ... 0) as unique eigenstate of eigenvalue 0, and all other eigenvalues are at distance at least sfKlL^I'^ if 
d < 2. H' can be regarded as a physical Hamiltonian. In Fig. |2] we represent H' for the case in which H is the 
frustration- free Hamiltonian of a system in one spatial dimension, and d—\. 

The construction in Ref. [T] leads indeed a family of frustration-free Hamiltonians in one spatial dimension. For a 
quantum circuit with Q gates, the spectral gap of iJ in Ref. [Ij can be made of order l/Q^ [H]. H' in Eq. (20) allows 
us to amplify the gap up to order XjQ^'^^l'^ ^ with an important gain if d = 2 and Q 3> 1. 



V. GAP: GENERAL HAMILTONIANS 



We show that spectral gap amplification is not possible for general Hamiltonians in the black-box model. Again, 
we consider instances whose unique ground states jV'a;) can be used to solve SEARCH as in Lemma [l] As before, 
the cost of an evolution with is the number of black boxes Ojj^ required to implement it. We obtain: 
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En. 



additional qubits 



FIG. 2: A physically local representation of H' . Each additional qubit 1, . . . , L has assigned a color and interacts with a 
particular term Hk in the Hamiltonian and the centered qubit. In a one dimensional construction, 11^ represents a nearest- 
neighbor interaction. Therefore, 11^ and Ilk+i have support overlapping on a single body that is represented by lines of color 
that match that of the additional qubit. 

Theorem 3 There are stoquastic ( sparse ) Hamiltonians for which spectral gap amplification is not possible. 

Proof 5 The result will follow again from the known lower bound on SEARCH. We consider the family {i?x}o<x<M-i 
with 



= -|a;)(a;| -c^ , 



(21) 



and X unknown. The constant c > will be specified later and A is a M x M symmetric matrix proportional to the 
adjacency matrix of a square lattice in five spatial dimensions. The matrix entries of A are Oyiy — 1 if {y,y') is an 
edge in the lattice and ay'y = otherwise. A is sparse and independent of the unknown x. In this case we assume a 
matrix oracle Wx for H^. This is stronger than Oh^ in the sense that Oh^ can be implemented with a single call to 
Wx. Further, by construction, Wx requires calling the search oracle of Eg. |Ip only once. 

We present some useful properties of Hx and refer to Ref. [TJI for details. There is a phase transition point at 
c — c* > where the spectral gap is 



A 



e 0(1/VA/) 



(22) 



Px = |(a;|V'a;)|' 

Ps^\{s\i:x)? 



0.37 



(23) 
(24) 



Then, the ground state \'4'x) has good overlap with both \x) and \s) as in Lemma\^ SEARCH can be solved with cost 
0{Vm) using the black boxes Oh^ or, similarly, the search oracles. If gap amplification for Hx were possible, this 
would yield a faster method to solve SEARCH. The proof follows by contradiction. 



VI. STOQUASTIC VS. STOCHASTIC MATRICES 

Because stoquastic Hamiltonians do not suffer from the so-called sign problem [37] . an important task is to devise 
classical probabilistic methods to solve the AST problem by sampling from the distribution as determined by the 
ground state. A well-known method for this purpose is quantum Monte-Carlo (QMC) '6]. Another method may be 
obtained by finding a diagonal similarity transformation that maps stoquastic Hamiltonians to stochastic matrices |13| . 
The fixed point of the resulting Markov chain is the desired distribution. Contrary to QMC, this second method 
requires knowledge of the ratio of amplitudes of the ground state for the mapping. Our result on no gap amplification 
for all stoquastic Hamiltonians implies that finding the transformation to stochastic matrices should be hard in the 
black-box model. Otherwise it would imply a method to amplify the gap by following the results on frustration-free 
Hamiltonians. 
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We can go further and give a family of stoquastic Hamiltonians for which classical probabilistic methods, that solve 
the AST by sampling from the ground state, has a large black-box implementation cost. For some QMC algorithms, 
our result implies that the spectral gap of the constructed stochastic matrix rapidly decreases with the problem size. 
Once more, our proof considers a family of Hamiltonians {-ffi;}o<rE<Af-i such that sampling from the ground state of 
Hx outputs the unknown x with large probability. 

Lemma 3 Consider the family of Hamiltonians 

Hx = -\x){x\+F/4, (25) 

where F G C^^ x C*^ is independent of x and \\F\\ < 1. Then, the lowest eigenvalue is bounded as Eq < —3/4, the 
spectral gap satisfies A > 1/2, and the ground state has large support on \x) , i.e. \{x\ipx)\'^ S ^(1)- 

Proof 6 (Sketched) All results can be simply proved by considering — |a;)(a;| + gF/A and using perturbation theory in 
the range [0, 1] for g. 

If F is stoquastic and irreducible, the ground state of is \ipx) = J2y \/^\y)j with iTy > 0. In this case we let 
Sx be a stochastic matrix whose fixed point is almost the distribution tt. We assume that a matrix oracle for Sx 
follows from a matrix oracle for Hx- The latter requires using the (classical) search oracle only once: for input y, the 
matrix elements of Hx in the yth row depend on whether x = y or not. In path integral QMC, for example, Sx can 
be constructed efficiently as follows. For sufficiently large /3 the probability TTy is (almost) proportional to 

{y\e~^"^\y). (26) 

Because A > 1/2, it suffices to choose /3 e ©(logM) to bound the approximation error. Next, because the expo- 
nentiation of Hx might be hard to compute, we write e'^^^" — nf=i^ ''^^i with rj G 0(1//?), and approximate 
Q-vH^ _ J „ _|_ (^(jj^y \Ye can relate nf=i(-'l ~ "H^x) with a "classical system" He (diagonal Hamiltonian) that 
acts on a larger Hilbert space of dimension . More precisely, we define 

with /r-i.r — 11 — jyi/a; I z*- ''-'). Notc that tr {e~^-^'') « tr (e^^"^). The computation of any diagonal entry 

of He requires evaluating p £ C'[polylog(A/)] matrix elements of Hx and thus the same number of search oracles. 
Sampling from the distribution determined by H^ requires choosing an appropriate Markov chain, like the one in the 
Metropolis algorithm. From known complexity bounds on the search problem, any choice of Markov chain will lead, 
in the black-box model, to a mixing time 



M 

polylog(M) 



(28) 



VII. CONCLUSIONS 



We introduced and studied the spectral gap amplification problem or GAP. The goal is to efficiently construct 
a Hamiltonian that has the ground state of a given Hamiltonian as eigenstate, but a bigger spectral gap. This 
problem is motivated by the results in Ref. j36j, where we showed quantum speedups of Monte-Carlo methods by 
giving a quadratic gap amplification for a specific family of Hamiltonians. The GAP has important applications in 
adiabatic quantum computing and other methods for adiabatic state transformations, where the implementation cost 
is dominated by the inverse gap: If gap amplification is possible, this results in a quantum speedup. 

We showed that a quadratic spectral gap amplification is indeed possible for frustration-free Hamiltonians. This 
generalizes the results in Ref. |36j and gives a more efficient construction to speedup classical Markov-chain based 
methods. In a suitable black-box model, we proved that our method provides the biggest amplification possible. The 
result for frustration-free Hamiltonians also provides more efficient quantum adiabatic methods to simulate quantum 
circuits and for the preparation of projected entangled pair states or PEPS. 

We presented a family of Hamiltonians for which spectral gap amplification is not possible under some assumptions. 
The Hamiltonians in this case can be sparse and stoquastic. They are constructed so that by preparing their ground 
states one can solve the unstructured search problem, and then use known lower complexity bounds to prove the 
impossibility result by contradiction. A corollary is that classical probabilistic methods that sample from the ground 
state distribution of some Hamiltonians may require a large implementation cost. 
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Appendix A: Proof of optimal amplification 



Each Hx is a symmetric matrix related to the adjacency matrix A of a sparse expander graph G = (M, d, A). M is 
the number of vertexes, d e 0(1) is the degree of G, and Xd is the second largest eigenvalue of A. WLOG, A < 1/2 
(e.g., a Ramanujan graph). Before specifying in a frustration-free form, we define by its matrix elements as 
follows: 



{y\H.\z) = { 



^ if {y, z} £ E and either y ^ x oi z = 

^ if {y, z} £ E and x ^ y ^ z ^ x 

1 li y = z ^ X 

otherwise. 



(Al) 



E is the set of edges of the graph G. has the following unique ground state of eigenvalue 0: 

1 , , 



71' 



v/2(M 



(A2) 



This state satisfies the requirements of the algorithm in Lemma [TJ \ipx) has constant overlap with \x) and the equal 
superposition state. Since x is unknown, the preparation of \iIjx) will allow us to solve SEARCH. 

We can bound the spectral gap A of Hj. in terms of M. We define Z^. = \ — — \tpx){''l^x\- Because \\Hx\\ < 1, 
the norm of is \\Zj:\\ = sup\^-^{4)\Zx\(j)) 1 - A. We can write |0) = I3\x) + ^/l - I3'^\x-^), with ^ > and \x-^) a 
unit state orthogonal to \x), to obtain 



mZxm < 2/3^1- l3^\{x\Zx\x^)\ + 
+ f3^\{x\Z,\x)\ + {l-fi^)\{x^\Z,\x^)\ 



An upper bound on the rhs of Eq. (A3) implies a lower bound on the gap. We write |a; ) — '^y^^lyi 
term on the rhs of Eq. ( A3 ) is 



y^x 



1 

2y/M- 1 



We use 



and 



to show 



{x\Hx\y) 



1 



2VM - 1 



< 



1 



2VM - 1 



V y 



y^x 



\2PVi^{x\Z,\x^)\ < \2f5Vl^\l 



(A3) 
The first 

(A4) 

(A5) 
(A6) 



(A7) 
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Further, from the diagonal matrix elements of , we have (M > 3) : 



P^\{x\Z,\x)\ < 



1 



1 



2 M ~1 



(A8) 



Finally, if Px — J2y^x \y){y\ projector onto the orthogonal subspace of \x), then the third term in the rhs of 

Eq. (IXsl is bounded by (1 - l3'^)\\PxZxPx\\- In addition. 



PxHxPx=Px{t-A/d)Px 



(A9) 



where A is the adjacency matrix of G(M, d, A). If \s) is the equal superposition state, A\s) = d\s) , and this is the only 
eigenstate of A with largest eigenvalue d. We can write 



PMx) 



M 



V^{M - 1) 



Px\s) , 



and consequently, 



< 



Px l^A/d 

A/d- 



M 



2(M - 1 
M 



-\s){^s\]P, 



2{M - 1) 

Because the second largest eigenvalue of A/d is A < 1/2, we obtain 

M 



s){s\ 



< 



Using Eqs. ( |A7[ ), ( |A8| , and (|Allj), we get 

((l)\Zx 



2{M - 1) - 2 



/?2 



^^^^^ M_l 



The maximum is found when /3 satisfies 



M - 1 (Af - 1)2 2(M - 1) 



Inserting Eq. (jAlSj) in Eq. ( |A12| , and since /3^1 - < 1/2, yields 

(01^x10) < 1 

< 1 

for M > 3. This implies 

A> ^ 

- 4(A/ - 1) 



1 1 

2{M - 1) ^ 2(Af - 1)2 
1 

A{M - 1) ■ 



(AlO) 



(All) 

(A12) 

(A13) 

(AM) 
(A15) 

(A16) 



We now focus on the specification of Hx (or a modified version Hx) as a frustration-free Hamiltonian. To this end, 
we first describe a particular choice for the projectors appearing in Hx- We define the function 



li y — X 



1 



(A17) 



We also define Eo as the set of ordered pairs (y, z) of edges in G, where each edge appears only once. That is, if 
{y, z} e E, then either (y, z) or (z, y) appears in E'q. For each edge (y, z) e E'o we define the (unnormahzed) state 



= c(y)|y) - c(z)|z) . 



(A18) 
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Then, 



H.= l</'(y,^))(</'(y..)K (A19) 

(y,z)eBo 



as the matrix entries coincide with those in Eq. (Al ). Because (0(^,2) = for all edges, is frustration free. 
For consistency with the definitions in previous sections, we define the projectors 



{4'{y,z)\4'{y.: 



n(y..) = 'T'^^ ,7"-^;' (A20) 



and the corresponding modified Hamiltonian 



"(y,.) . (A21) 



\ipx) is still the unique ground state of . Since 



{^\Hx\ip)= J2 i^\^{yM> i^\hy,-))i^{yM = ' (A22) 

{y:z)eE„ {y,z)eE„ 



the spectral gap of Hx is also bounded from below as A > 



4(M-1) ■ 



The number of terms in Eq. ( A21 1 grows with the number of edges in the graph, dM/2. It is convenient to reduce 
the number of projectors so that L £ 0{1), and L playing no significant role in the scaling of the amplified gap. We 
then assume a given edge coloring of the graph G with edge chromatic number x of order d |41j [i.e., x G 
Let ci, . . . be the different colors. All the projectors ^{y.z) belonging to one of the colors are, by construction, 
orthogonal to each other as they don't share a vertex. For each /c G {1, . . . , x}, we define the projectors 

Y , (A23) 

Then, i?, =ELinfc. ^ 

The specification of in this proof is as follows. For each k we assume the existence of a reversible version of 
the matrix oracle W that, on input {k,y), it outputs (z, (z|nfe|j/)), where z is the only vertex that shares an edge of 
coloring Ck with y. Note that all the matrix elements (zjllfel?;) are equal unless y = x or z = x. The implementation 
of W requires then deciding whether its input y is the marked vertex x or not. It implies that the reversible version 
of W can be implemented using a single call to the search oracle of Eq. (jT|. Further, the black box , which 
implements exp{— jn^i} on ancilliary input can be constructed using 0(1) oracles W. This follows from the 

simple observation that exp{— illfct} induces rotations in the two-dimensional subspaces spanned by {\y), \z)}, where 
z is the only vertex that shares an edge of Cfc with y. The implementation of a reversible requires then using the 

7 7 

search oracle at most two times (y = x or z = x). 

Let H'^ be the Hamiltonian that has | V'^)|o) as a non-degenerate eigenstate and spectral gap A'. SEARCH can then 
be solved using the technique in Sec. O which requires evolving with H'^ for time 0(1/A'). By assumption, it implies 
that SEARCH can be solved using 0(1/A') black boxes or, similarly, 0{1/A') search oracles. Therefore, the 

known complexity bound on SEARCH implies that A' < XjyfM e C'(-\/A). In this black-box model, our construction 
in Sec. |III| gives the biggest gap amplification possible for frustration-free Hamiltonians. 



Appendix B: Spectral properties of G 



The first property to notice is that if |^o) is in the null space of i?, then H^j'^o) = 0. Therefore G'|'0o)|O) = and 
|V'o)|0) is an eige nstate in the null space of G. 

As in Sec. HI we label the eigenvalues of H by \j. Assume now Aj ^ 0, and consider the action of G on the state 
1^,0) 



(Bl) 
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Notice that 

(V,-0|G|^jO) = ^(V',fc|nfc|V',0) = (B2) 

k 

and that 

\mm' = E(^^-^infclV'.fc) = A, . (B3) 



We denote by 



1^.) = ^G|V,0) (B4) 
V -^J 



the corresponding normahzed state. Next note that 

G\±j) - ^ ^^n^i^.o) - v^lV.O) . (B5) 



Therefore, the Hamihonian G is invariant in the subspace Vj = {j-f/'jO), |-Lj)}. Define Gj = G|y^. the projection of G 
in this invariant subspace. In the basis {iT/jjO), |-Lj)} we can write 



A, 
with eigenvalues ±-y/Aj. 

Finahy, we note that because the form a complete basis, G acts trivially on any other state orthogonal to 
(BjVj- While the dimension of the null space of G is larger than that of the null space of H, this has no effects when 
solving the AST problem: transitions between |V'j)|0) and any other orthogonal state in the null space are forbidden. 
We can also use the trick in Sec. |III| to reduce such a degeneracy, if necessary. 
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